Mode coupling evolution in arbitrary inflationary 
backgrounds 



Francis Bernardeau 

E-mail: francis.bernardeauOcea.fr 

CEA, Institut de Physique Theorique, 91191 Gif-sur-Yvette cedex, France 
CNRS, URA-2306, 91191 Gif-sur-Yvette cedex, France 

PACS numbers: 98.80.-k, 98.65.-r, 98.80.Bp, 98.80.Cq, 98.80.Es 



Abstract. The evolution of high order correlation functions of a test scalar 
field in arbitrary inflationary backgrounds is computed. Whenever possible, exact 
results are derived from quantum field theory calculations. Taking advantage of 
the fact that such calculations can be mapped, for super-horizon scales, into those 
of a classical system, we express the expected correlation functions in terms of 
classical quantities, power spectra, Green functions, that can be easily computed 
in the long-wavelength limit. Explicit results are presented that extend those 
already known for a de Sitter background. In particular the expressions of the 
late time amplitude of bispectrum and trispectrum, as well as the whole high-order 
correlation structure, are given in terms of the expansion factor behavior. 

When compared to the case of a de Sitter background, power law inflation 
and chaotic inflation induced by a massive field are found to induce high order 
correlation functions the amplitudes of which are amplified by almost one order 
of magnitude. These results indicate that the dependence of the related non- 
Gaussian parameters - such as /nl - on the wave-modes is at percent level. 



1. Introduction 

It is now clearly established that standard single held inhation cannot produce 
signiheant non-Gaussianities (NG) during or immediately after the inflationary phase. 
This has been explicitly shown by Maldacena in Ref. [1] where it clearly appears that 
standard single field inflation leads to no or very little primordial non-Gaussianities. 

Multiple-field inflation has now long been recognized as a possible mechanism for 
the generation of primordial metric NG fluctuations (see recent review papers in Refs 
[2, 3, 4, 5]). The exploration of the various types of coupling terms that appear in the 
action then leads to distinguish gravity from non-gravity mediated couplings (see [5] 
for recent review on the origin of this distinction) . This decomposition comes from the 
various behaviors of the terms that are present in the third order action. In one-field 
inflation, because the field fluctuations and the metric fluctuation are locked together, 
only the former can be found. In multiple-field inflation however this is not necessarily 
so. In particular isocurvature degrees of freedom open the possibility of having a 
richer phenomenology. And whereas gravity mediated couplings are ubiquitous but 
induce only modest effects [6, 7], non-gravity mediated couplings can be very efficient, 
although nothing ensures that they are generically at play (this is at the heart of Refs. 
[8, 9, 10]) . 
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To make such perturbations play a role, one indeed needs a mechanism to transfer 
isocurvature modes into adiabatic fluctuations. For instance the curvaton model is 
based on the survival of (massive) isocurvature modes until late after the end of 
inflation that can alter the subsequent expansion history of the universe [11]. This is 
a particular case of modulated inflation [12, 13, 14]. Other mechanisms assume that 
isocurvature modes can change the end-point of inflation or alter the (p)-reheating 
sequences (as in [15]). Such mechanisms can also happen in the context of hybrid 
inflation. This latter situation is in particular advocated in Refs. [16, 17, 10, 9] 
where isocurvature modes are shown to be able to induce large NGs in the metric 
fluctuations. 

In this paper we assume that such a mechanism is at play and focus our analysis 
on the intrinsic statistical properties of isocurvature modes, that is of modes that do 
not participate in the metric fluctuations at the time we are interested in. Even though 
this is a simple setting, exact results are still difficult to obtain since they depend on 
the actual dynamics of the expansion of the universe. Some exact results have been 
obtained in case of a de Sitter background (in Ref. [18] bispectrum of a field with a 
cubic potential is derived; in Ref. [1] more complicated types of cubic interaction are 
considered and in Ref. [19] the four-point function is derived for quartic interactions) 
that has been the favorite playground for such calculations. The evolved bispectra, 
trispectra, have in particular been found to grow like the number of cfolds after horizon 
crossing. This is an important result that lays the foundation for the computation 
of the mode dependence of the coupling parameters (which are usually expressed in 
terms of /nl parameter as introduced in [20] ) . This possibility is now being explored 
(as in Ref. [21] for instance for some classes of model.) However, though inflationary 
models lead to background evolution that are generically close to a de Sitter, this is 
only an approximation, and which at best can be valid for only a limited period of 
time. 

In this article we aim at being as much exhaustive as possible in presenting exact 
results on the evolution of the statistical properties of a test scalar field, with non- 
vanishing self interaction terms, in arbitrary inflationary background. We will make 
use of the relation between quantum and classical evolutions. However we will not use 
explicitly the SN formalism (as introduced in [22, 23] and see Ref. [24] in the context 
of nonlinear expansions) as it does not necessarily give controlled approximations but 
try to be as much precise as possible to draw the line between results that are of 
sub-Hubble origin and those that can be accounted for in a super-Hubblc classical 
calculation. In the "tree-theorem" in Ref. [25] Weinberg demonstrated that there 
exists a classical system that exhibits the same tree order correlations properties as 
the quantum system. Although we will not exactly use this solution, it will eventually 
lead us to a complete description of the super-Hubble tree order correlation functions 
of the field. 

The paper is divided as follows. In section 2 we recall the method and the known 
results in case of a de Sitter background and describe in some details the late time 
super-Hubble limit of bispectra and trispectra. In the following section we make 
explicit the connection between calculations of correlation functions in a quantum 
context and classical calculations of stochastic field evolution. It shows in particular 
that the late time evolution of correlation functions can be computed in a classical 
context with the introduction of the Green function of the classical evolution of the 
free fields. This observation is exploited to obtain new results for the computation 
of the late time behavior of correlation functions in arbitrarily backgrounds. These 
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results are presented in the section 4. In section 5 the resulting correlation functions 
are presented in a systematic way with the description of the tree structure of the field 
correlation functions at leading order in perturbation calculations. 



2. Mode coupling computation 

2.1. General framework 

We are interested here in the evolution equations governing a test scalar field x( x , v) 
in the homogeneous background of an inflationary universe. Throughout the paper 
we assume that the spatial curvature vanishes so that the background in which the 
field evolves can be described with the following metric, 

ds 2 = a 2 (r;)(-d77 2 +da; 2 ) (2.1) 

and a{q) determines the background expansion. The time dependence of 0(77) is 
assumed to be given by a sector of the theory which is independent of \- We 
neglect in particular here gravity induced loop effects such as the backreaction of 
the inhomogeneities in \ on the background evolution. The \ field therefore behaves 
as a minimally coupled test field. We further assume that the interaction operator of 
X, V [x], is a function of \ only. This is again for convenience since the method we 
introduce could be a priori extended to more more complicated operators!. The aim 
of the following calculation is then to explore the statistical properties of the x field 
or more precisely of its Fourier modes Xk(^) 7 

Xk(v)= y d 3 x exp(-ik.x)x(x,J?), (2.2) 

that is compute the late time power spectrum of the latter, V x (k), defined as 

(XkXk<) = (2^) 3 Mk + k')P x (fc) (2.3) 

where is the 3-dimensional Dirac delta function and its higher order correlation 
functions such as the bispectrum, _B x (ki, k 2 , ks), defined as 

(XkiXk.Xks) = (27r) 3 <5 D (k 1 +k 2 +k 3 )B x (k 1 ,k 2 ,k 3 ). (2.4) 

For convenience we also define the reduced high order correlation functions such as 
Q 3 ({h}) or Q 4 ({fci}) generally as 

<Xk 1 ...Xk p ) c = (27r) 3 <5 D (^k l ) Q P ({fcJ) EIl P x(%) (2-5) 

and our final results will be expressed in terms of Q3 or Q4 functions for simple 
potentials. 

For the actual resolution of the equations, it is convenient to introduce the reduced 
field u(x, 77) as 

u(x,r?) = a(i)x(x,77). (2.6) 

The resulting motion equation is then, 

a" dV 

«"(x, 77) u(x, r?) - A 7i(x77) = -a 3 — [u(x, r?)/a] (2.7) 

a dx 



I As this is the case for the cubic part of the action for the inflaton field for instance. 
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The right-hand side of this equation contains nonlinear terms in u(x, 77) and this 
evolution equation cannot then be solved explicitly. One can then define the free field 
u(°) which is the field that satisfies the motion equation when only the quadratic terms 
in the Hamiltonian are kept. 

In the following we put an " over u, or on other fields, when we insist on the fact 
that it can be viewed as quantum operators. It is then assumed that the field u 1 - ' can 
be decomposed into creation and annihilation operators, ak and such that 

fl(0) (x ' v) = j (S* exp(ik - x) ^ 0) {r]) (2 - 8) 

with 

a k %)-4%w+4 0) W- k (2.9) 

where we have 

[a k ,al] = ft(2^) 3 Mk-k'). (2.10) 

In Perturbation Theory one makes the assumption that the field u(x, rf) can be 
expanded in terms of those same operators. The time dependence of the function 
u k°^ (v) i s set by the time dependence of the background. In the following we will be 
interested in cases where the mass of the free field vanishes (or is assumed to be small 
enough compared to the Hubble constant) so that the motion equation is, 

uf'irj) ^4%) + fc 2 4%) = 0- (2-11) 

In principle this equation can be solved once a{n) is known. 

In the following, we will then assume a simple form for the potential. 

V( X ) = ^X P (2.12) 

where p is an integer, p > 3, which in practice will be cither 3 or 4 corresponding to 
respectively a cubic or a quartic coupling. 



2.2. Conditions for large isocurvature fluctuations 

It is probably worth investigating in some details the conditions (on A) for x to develop 
fluctuations in the first place. A priori such a potential is bound to induce an effective 
mass square to the field % of about \H V ~ 2 so that one should have A < H 4 ~ p in 
order to fluctuations in the x-field to develop. For a cubic coupling that implies that 
A < H and A < 1 for a quartic coupling. The derivation of this constraint can be made 
more quantitative with the help of a classical stochastic approach (introduced in [26] ) 
which can be used to infer the probability distribution function of the value of x a t 
horizon crossing. The Fokker equation for the time dependent probability distribution 
function of x, V(t, x) can be derived from the evolution equation of x- It is given by, 

&P = H 3 d*V 1 d (dV{ X ) \ 

dt 8tt 2 d X 2 3H d X \ d X J ' ' 

where t is the physical time. This equation has a non trivial late time, and time 
independent, solution given by (for p = 4), 

7T 2 A X 4 - 



rix) = wwm It l exp 



9H 4 



(2.14) 



Mode coupling evolution in arbitrary inflationary backgrounds 



5 



It shows that the typical excursion values for x are 01 the order of A _1 / 4 i/ for a quartic 
coupling. It corresponds to excursion values of a field of mass of the order of \ X I 4 H. 
Modes of a given wavelength, whose amplitude if the field is effectively massless, are 
therefore not suppressed if A < 1. 

This result can be further comprehended from the expression of loop corrections 
to the two-point correlation of self coupled a test field. In [27] it is shown in particular 
that a complex scalar field with a quartic coupling imbedded in a supersymmetric 
multiplet developed an effective mass square of the order of AH 2 times a logarithmic 
divergent term (which in the context of inflation is reasonably finite.) This is, if one 
needs to be fully specific, the context that we have in mind in the investigations we 
present here. Note that this condition ensures that the isocurvature fluctuations are 
of the order of H at the time of horizon crossing. It does not prejudge the conditions 
for the survival of the isocurvature perturbations till the end of inflation (where it is 
assumed that some sort of mode transfer can take place.). For de Sitter background 
perturbation theory results suggest that the condition A < H 4 ~ p should be replaced 
by XN e < H 4 ~ p where N e is the number of efolds between horizon crossing and the 
end of inflation. We will see at the very end of the paper how this condition is to be 
(slightly) changed for other backgrounds. 

2. 3. The "in-in" formalism 

One unambiguous way to derive the statistical properties of the field x is to use 
the "in-in" formalism from Schwinger and Keldysh (see Refs. [28, 29] and [30] for a 
presentation and use in a cosmological context). We briefly present its implementation 
in this context here. 

The first step is to decompose the Hamilonian of the system in a classical part, 
the one responsible for the background evolution, and a quantum part that will govern 
the evolution of the quantum fluctuations. The latter contribution can be further split 
in two parts, 

H Q =Hf+H^\ (2.15) 

one which is quadratic in the fields, Hq \ and which governs the evolution of the free 
fields, and an interaction term, , which is supposed to be small so that the actual 
evolution of the fields remain close to the free field evolution. This is obviously a basic 
hypothesis for doing perturbation calculations. This decomposition allows to define 
the operator, U^\ri, rjo), as 

Aw« ( V , Vo ) = -iff « [ V , x (0) (v)} (ry, no) (2.16) 

where the interaction Hamiltonian is written in terms of the free fields. The solution 
of (2.16) formally reads, 

«< J > fa, m ) = T jcxp (-i j" [ v , X (0) (r?)]dr?) j . (2.17) 

In this expression the operator T is the so-called "T- product", meaning that the 
operators that appear in its arguments should be ordered in time. Specifically the 
terms that appear in the argument of the exponential should each be ordered in time, 
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e.g. 



J no 
rn pn 

- d m d m H^( m )H^( m ) + ... (2.18) 

J Vo J no 

Then using this operator one obtains for any operator Q(r]), that is any function of 
the field x(x, 77), 

Q{rf) = (U^y'^Vo) Q (0) (l)W (J) (l,%). (2.19) 
which again can be written, 



Texp (\ f' H {I) dr]' 

\ J Tin 



Texp -i 



H^dr,' 



•(2 



This relation is the keystone of perturbation calculations for observable in an 
expanding universe. It implies for instance that 

(0|Q|0) = (0|Q (0) |0) 



-i P dm (0| [Q (0) , H^( m ) 

J rin 



dm 



no 
+ ... 



d m (o| [[q (0 \ H^(m)} ,H^(m)} |o) 



(2.21) 

Each line of the preceding expression corresponds to a given order in perturbation 
theory. It is easy to see that higher order terms will involve expression with more and 
more entangled commutators. § 

Let us illustrate this property for observables involving simple field products and 
for an interaction hamiltonian that takes the fom, 



H^\r})dr] = a^^dr] J d 3 x V u (0) (x, v)/a(v) 



(2.22) 



where V is of the form of Eq. (2.12). In the following, we will not in particular consider 
terms that involve higher order derivatives of the field (although the existence of such 
terms cannot be excluded) and in practice we will consider terms that are cubic or 
quartic in the field. 

As can be easily seen from the terms involved in (2.21), the whole result can be 
written in terms of the unequal time correlation function of the free fields defined as, 

(0|u k (77)u k ,(7/)|0) = (2^) 3 <5 D (k + k') G k (mv')- (2-23) 

In can be expressed in terms of the function u^(rj) that are solution of (2.11). And 
we have 



G k (mv') = 4%K u; V)- 



(2.24) 



§ There are alternative formulations of these expressions in terms of time dependent propagators, 
see for instance [31, 32]. 
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In the following we explicitly give the expression of such correlation functions in 
case of a quartic potential. From the relation (2.21) we have, 



(Xki • • • Xk 4 )c — 



1 



a 4 ( V ) 



(u kl . . . Uk 4 )c 



cV(0| [«£> ( V )...u^( V ), Hj{rf)l |0) 

/' d?? '(2^ dk --- dk4 ' :X N; 
.4%),4?(V)..-4?(V)1|0> (2.25) 



iA 



4! a 4 (r?) 



x (Olfu^W 



It finally leads to, 



(Xki ■ • • Xk 4 )c — 



iA 



(2^) 3 <5 D hTk, 



x / d v '[G kl (v,v')---GkAv,v')-cc.} (2.26) 

where c.c. stands for the complex conjugate of the previous term. Note that from the 
definition of Gk{f],rf) in Eq. (2.23) it is clear that Gk(r]' ' ,r}) is the complex conjugate 
of Gk{f],rf). In practice the computation of such cumulants relies therefore on our 
knowledge of the explicit form of the G functions. In the following we present two 
specific cases where the calculations can be completed analytically. These are the 
cases of de Sitter background and some models of power law inflation. 

2.4- Three- and Four-point function in de Sitter background 

-2/rf and the motion equation of 



For a de Sitter background we have a" (rj) / a(r]) 
Eq. (2.11) can be solved exactly. We then have 

„(%) = J= (l- * ' 

It is interesting to note that this solution can be rewritten, 



ikt] 



4%) = 



2k rj 



where the operator Ok is defined as 

o k = 1 --^. 

k k dk 

As a result, in this case, the expression (2.26) can be rewritten, 

/ 4 



(2.27) 
(2.28) 
(2.29) 



(X^ • • -Xk 4 )c 



iA 



I6ir 4 r) 4 a 4 (ri) 
O kl ...O k4 e-^A (o kl ...O k 



(2tt) 3 <5 d PTk, 



\i=l 

" dr/ 
n 



14 



(2.30) 



where we introduce the following combination of the four wave modes, 

7Tl = y^fci, 7T 2 = ^2 ki kj , 

i i<j 

7T3 = ^2 ki kj kk , 7T4 = ^2 ki kj kk h . 

i<j<k i<j<k<l 



(2.31) 
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The integral appearing in Rcf. (2.30) reads 



/ 



" dry' 



1 



7 i Ei (i?77ri)7rj + 



e i^ l7r 2 i e i^ l7ri £ i^ 1 



(2.32) 



r;"* 6" -"-'■■ ± ' ' 6r? 6r? 2 3r? 3 

from which the explicit expression of (xki ■ • • Xk 4 )c can easily be derived (it is given 
in Rcf. [19] in terms of the combinations 7r,). 

In the the super-Hubble limit (i.e. kit] <C 1 for alH = 1 . . . 4) the expression of 
Q4 as defined in Eq. (2.5) can be explicitly computed giving, 

Q 4 ({M) - ^2 + + lo § (-^E k *)} ■ ( 2 - 33 ) 

In this expression, terms of order fcj?y have been neglected. The £ function is an 
homogeneous function of the wave vectors which reads, 

' ' 2 7Tl 2 7T2 + 7Tl 7T3 — 3 7T4 



C({fcJ) 



-7T1 



(2.34) 



7Tl (7Tl 3 - 3 7Tl 7T 2 + 3 7T 3 ) 

It depends on the wave vector ratios but in a rather weak way. Its minimum value, 
—51/16, is reached for a symmetric configuration of the 4 wave modes. Its maximum, 
—2, is reached when 2 of the wave mode values vanish. 

Obviously similar results can be obtained for the 3-point function in case of a 
cubic potential. In this case the quantity to compute is 



(Xk! •••Xk 3 ) 



O kl ... O k3 e 



iA 



87r 3 ?7 3 a 3 (?7) 



(2^) 3 <5 D k * 



-17Ti>j 



o. 



i=l,3 
v ,4 e 



(2.35) 



Q 3 = 



3H 2 



Interestingly, to a factor 2n, one recovers the same expression as for the four-point 
function when 7T4 is put to zero so that Q3> once again defined from Eq. (2.5), takes 
the form, 

r 7+c 3 ({fc,})+iog(-'?E^)] ( 2 - 36 ) 

in the super-horizon limit where (3 is, 

&({**}) = C(fci, to, fc3,*4->0). (2-37) 

The values of C 3 ({^}) ran ge from —8/3 to —2. This expression gives the ki dependence 
of the usual parameter for the x field - which is nothing but Q3/2) - as, 



A 



6H 2 
kik 2 k 3 



N e (k 3 ) - log 



fc 3 



■7 



(fci + fc 2 + fe 3 ) (fc 2 + k\ 



j_ Z~3 j_ Z~3 



(2.38) 



were N e (kz) is the number of efolds between the horizon crossing of k 3 and current 
time. A similar result can be derived for = Q4/6. They both scale like XN e /H 2 
and therefore give an order of unity parameters for AiV e w _ff 2 which would then be an 
effect comparable to the gravity induced couplings. Those results are therefore valid, 
and of interest, for H 2 < \N e < H for a cubic coupling and H 2 < \N e < 1 for a 
quartic coupling. 
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Figure 1. Dependence of on fcj as a function of k\ and &2 for a fixed value of 
&3 and for X/H 2 = 1. This function is computed assuming the numbers of efolds 
since k$ crossed the horizon is N e = 50. 



The resulting fc— dependence of is shown on Fig. 1 for a N e {k^,) — 50 
and X/H 2 = 1. What this figure shows is the remaining dependence of the non- 
linear coefficient on configurations. Clearly the main feature is that /nl is almost 
independent of the configuration and of scale. To be more precise it is straightforward 
to see that, for a given configuration (e.g. given ratios of fcj), 

dl °S^NL = 1 f2 39) 

dlogfc t N e (k t ) [ ' 

where k t = ki + fc2 + k^ and A^e is the number of efolds since kt crossing. We are thus 
left with a 2 percent effects. On Fig. f one can see however that the dependence of 
/^j L on the geometry can be as large as 7.5 percent for wave mode ratio of fO. These 
dependences are entirely encoded in the form of Eq. (2.38). In the following we will 
see how they depend on the background dynamics. 

2.5. A case of power-law inflation 

There is another simple (yet unrealistic) model of inflation which can be solved 
explicitly and that corresponds to a case of power-law inflation. In the latter case 
the expansion factor is assumed to behave like t v with v > 1 where t is the physical 
time or equivalently a(rj) ~ I n such cases the evolution equation for u ® (77) 

can be solved in 



4%) = ^4 2) (- fc ^ (2 - 40) 
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where is the function of Hankel of the second kind with 

H=\ + ^-r- (2.41) 
2 v — 1 

When \i is a half integer it is possible to explicitly compute the expression of the 
four-point function of \- Besides the dc Sitter case (that corresponds to \i = 3/2), the 
case 1 1 v = 2 leads to [i = 5/2 and for which we have 

ui ° )(v) = v^f Qke ~ ikn {2A2) 

with 

3 3 d d 2 

Qk = ^-kOk + o^- (2 - 43) 

In this case the late time power spectrum reads 

P (fc) = ? I = ?MM (2 44) 

where Hq is the value of H at time of horizon crossing for the mode fco- The formal 
expression of high-order terms can also be computed. For a quartic potential it leads 
to the following form, 

(Xki ■ ■ ■ Xk 4 )c = - -z ^ 47 ^ (27r) 3 fe [ V k 4 

167T 4 ?? 4 a 4 (?7) ^.f^ 4 

x | ^Q fcl . . . Q fe4 e- iT1 ^ ^Q fcl ■■■Qk i f^ e^ 1 "') - c.c.j . (2.45) 

The super-Hubble limit of this expression can again be computed at leading orders in 
super-horizon. It leads to 

A a 2 (r])r] 2 



Qa ({h}) 



14 

1+-^ — 15EfcF log( -^^ 



. (2.46) 



The growth of the amplitude of the reduced correlation function is then extremely 
rapid since it grows like a(rj), that is like exp(iV e ) where N e is the number of efolds. 
The corrective term is only in N e . This is clearly at variance with the result in case 
of the de Sitter background case. That implies in particular that the k dependence of 
the and parameters is completely suppressed, by many orders of magnitude. 

In the following we aim at exploring in a systematic way the dependence of the 
late time behavior of the function Q p on the background dynamics. It will not be 
possible however to do so from exact quantum calculations with the "in-in" formalism 
and we rather take advantage of the classical re-formulation of the mode coupling 
evolution for super-horizon scales. 



3. Quantum to classical descriptions 

The classical system we introduce now is not exactly formulated as in Ref. [25] . It 
will nonetheless exhibit the same super-Hubble properties. 

|| There are actually an infinite number of cases that could be solved following this procedure but 
they would correspond to smaller values of v. 
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3.1. A classical stochastic counterpart system 

So let us now introduce a classical system that consists in a stochastic real field, 
u(x, if), the time evolution of which is given by Eq. (2.7). As a classical field, the 
amplitude of u(x, rf) is determined by its initial conditions that we take corresponding 
to those of the free quantum field u(°)(x, if). 
Remark that if we define 

§k = (ak + a^ k ) and d k = (a k - a^^j (3.1) 
then the relation (2.9) can be rewritten 

* (0) M = / (^ eik ' X [Rc(4° ) (r / ))s k + Im(4 0) (r ? ))d k ] . (3.2) 

We can then note that 

(0|s k s k <|0) = (0|d k d k ,|0) =ft(2^) 3 ( 5 D (k + k'). (3.3) 

For the classical model, the free field u^ ^ obeys Eq. (3.2) when § and d are "degraded" 
into time independent classical random variables that form each a Gaussian field. More 
precisely we introduce the random variables s and d as time independent Gaussian 
random variables whose ensemble averages are given by, 

(s k s k /) = (d k d k -) = h (2^) 3 fc(k + k'), (3.4) 

so that 

u (0) (M) - J [Re(4 0) (r / ))s k + Im(4 0) (r ? ))d k ] (3.5) 

is a classical free field the evolution of which is fully compatible with (2.11) (this 
is clearly the case since both Re(u^(r/)) and \m.{u^\rf)) are solutions of the 
linearized field equations). As a result the classical modes, u^\r]) = Re(u^\rj)) s k + 
lm(u^\ri)) d k , obey the following properties, 

(uf'Wul!^')) = (2nfS D (k + k')P k (r 1 , V ') (3.6) 

with 

P k ( V ,r,') - Rc(4 O) ( ?? ))Re(4 O) ( ?? 0) +Im(4° ) (7 7 ))Im(4 0) (r / ')) 

= Rc(4 0) ( ?? )4 0) *(r / ')) (3.7) 

which is nothing but the real part of Gk(i], r}'). 

The connection between the two systems is actually more profound and extends 
to non-linear evolution properties. It can be easily checked that the imaginary part of 
Gk{f],rf) is - to a factor 2 - the Green function of the classical system. The classical 
Green function G c (r],r]') can be written 

G^(v,v') = e(v-v')gM) (3.8) 

where 8 is the Heavyside function and guiViV') can De constructed from any two 
solutions of Eq. (3.2). It can easily be shown that, given the normalization constraints 
on (77) , we have 

G k ( V ,r}') = P k ( V ,r)') - '-g^n'). (3.9) 

It gives some direct insights into the transition of the quantum evolution into the 
classical regime. Indeed, at super-horizon scales, the phase of u k °\ri) tends to be fixed 
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which implies that the imaginary part of the G k {j],r)') functions decays with respect 
to its real part. As we see here this behavior is closely related to the classical evolution 
of the field for super-Hubble scales since G k (r], 77') then tends to be simply the classical 
unequal time field correlation function. 

This observation has also some consequences on the expression of the r.h.s of Eq. 
(2.26). Indeed it can be noted that as one takes the imaginary part of a product of 
Gk functions, the result necessarily involves at least one g k factor. In other words, the 
expression (2.26) can be rewritten, 

(Xki • • • Xk 4 )c = 




a*(r?) 

x f dr?' [P kl (77, 77') . . . P k3 (77, r,')g k4 (77', 77) 

-^ft 1 (?7,?7').9fc 2 (»7,77')---5fe4('?,»7')+ s y m - • (3-10) 

Naturally the first term appearing in the bracket is to be the leading contribution for 
super-horizon scales and it will turn out with a simple classical interpretation. 



3.2. Classical evolution of the high-order correlation functions 

We focus here on the mode coupling effects induced by simple forms for the potential, 
that is (iikj . . . 7ik 3 ) in case of a cubic coupling and (u^ . . . u^ 4 ) c in case of a quartic 
coupling. In order to be able to apply a perturbation theory approach we have to 
assume that can be expanded in terms of the coupling constant, 

«k(77)=4 0) (77)+4 1) (r7) + ... (3.11) 

where 7t<") is of order n in the coupling constant. The first order expression of u can 
easily by computed from the classical motion equation, e.g. we have, 

£ ) \v)-^u£\v) + #$\ri) = -Av) / d 3 x^e- ik - (3.12) 

where ^ is the derivative of the potential with respect to \ an d taken at value 
(x, 77) / 0.(77) . The solution of this motion equation can be written with the help of 
the classical Green function, g k (r],ri') as > 

v£\ V ) = - f a 3 (v')9k(v,v') I d 3 xcxp(-ik.x) ^. (3.13) 

In the following we treat the case of a quartic^f coupling V(x) — |fX 4 - In this case, 
the double Fourier transforms lead to, 

~ 3!(2^ f 9kM) J d3k i'-- d3k 3 

x4V)---4%0fc(k-£lq). (3.14) 
^ the case of a cubic coupling can be treated very similarly. 
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As a result, the tree order contribution of the fourth order cumulant reads, 



(Xk! • • • Xk 4 )c = - ^4(^0 (^E k * J 

x/ dr? Pk 1 (v,v')Pk 2 (v,v') p k 3 {v,v') 9k 4 (v,v') + sym. (3.15) 




This expression reproduces Eq. (3.10) when the last contribution has been dropped. 
At variance with the classical system proposed in Ref. [25], the evolutions of the two 
systems are here not identical, even at tree order. However, provided the integral 
appearing in this expression is dominated by super-Hubble contributions - and this 
is indeed the case for the two examples we know how to compute - so that gki/lirf) 
is much smaller than Pk(f],n'), the two contributions should be nearly equal. This 
is this observation that allows us to sort out the leading contributing terms to such 
correlation functions. 

4. Dependence on the background dynamics 

4-1. General expressions 

The practical computation of the leading behavior of Eq. (3.15) can then be further 
simplified by observing that one can use the classical Green function of the system in 
its limit + k — >■ 0. Furthermore, for superhorizon scales, the dependence of Pk{r),rf) 
on T}' is that of a(r/). The expression (3.15) is then simply proportional to, 



It is interesting to note that a calculation with a cubic potential would actually have 
led to the very same expression for Q^{r], {h}). In the following we will then note 
^eff.(^) the common values of Q p {rf) in the super-Hubble limit. 

We now turn to the explicit computation of goiVtV')- As we are considering 
massless fields, it can easily be derived from the behavior of the expansion factor. 
One solution of the field evolution is indeed given by, vi(n) — a(r)), the other solution 
is given by V2(r]) = a(rf) p A-q' / a 2 {rj'), so that the classical Green function then reads 
QoiVirf^iv ~ V r ) m the fc — »■ limit with, 




(4.1) 



which is nothing but the super-Hubble limit of Q4, 



(4.2) 




(4.3) 



The expression (4.1) can then be formally written, 




(4.4) 



in terms of the expansion factor. This is one important result of this paper. 



+ It is to be noted that the sub-dominant terms cannot be accounted for in a perturbative expansion 
of g with respect to kr\. For instance next-to-leading order in de Sitter case can only be derived from 
a full use of Eq. (3.10). 
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In the following we are going to explore the dependence of this quantity with the 
background dynamics. Let us start with the de Sitter case for which we have, 

gf(v,v') = ^(v 3 -v' 3 ) (4.5) 

which is indeed what the late time expansion of the imaginary part of G ds gives. We 
can now plug in this behavior into (4.1). In this case we then have, 

During the super-horizon evolution 77 gets close to —0 and 77/770 gets very small. The 
number of efolds after time 770 is nothing but N e = — log(r]/r]o). It corresponds, to 
a finite correction, to the number of efolds after horizon crossing, e.g. the validity 
domain of the form (4.3) of gk(rj,r] r ). In case of a de Sitter background the vertex 
value is then i/f (r?) = XN e /3H 2 . 



4-2. Power-law inflation 

We can now return to the power law case but without any restriction on the index. So 
we assume that the expansion factor behaves like t v with v > 1 where t is the physical 
time. In this case, we have H = vjt and 

ft Jj/ iV 4-IV 

3o Pl (M') = ft* J ^ = -^~ x {t 1 ^ t'^) (4.7) 



which implies that 

4-liv 



o df ^ (t^ - t'^) 

x (t 2 -t 2 t 2 -tl v+1 /t 3v - 



(4.8) 



(4.9) 



Zv- 1 V 2 3^ + 1 

For v = 2 one recovers the late time behavior found in the exact calculation. In 
general this expression can be reformulated in terms of N e and Ho, the value of H at 
horizon crossing time to. We have t/t = exp(7V e /^) so that, 

Xv 2 



1 ( 2N e /u _{] _ 1 ( e 2N e /v _ (l-3 V )N e /u\ 

2 V J 3v+ 1 V J 



(4.10) 



so that 



2 V J -Zv + 

If one further assumes that v/N e is large enough, then we have 

iVe = l0g^=l0gffy«^ (4.11) 

a(to) \toJ t 

^)«^(*-^»^, (4-12) 

which reproduces the result of the de Sitter case. Note however that this convergence 
is very poor. A simple expansion of (4.10) gives 

pi / \ XN e / N e 2N 2 \ ,„ 

^ W ^3^( 1+ 3i7 + 3# + -j- (413) 
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If we want to put realistic numbers for the index v, we can use the constraint 
provided by the value of the spectral index n s of the CMB anisotropics (as summarized 
in Ref. [33]). It leads to v « 40. Such a value invalidates the previous expansion. On 
the left panel Fig. 2 we then compare the de Sitter case (lower curve) to the power 
law inflation case for v varying from 35 to 50 (that corresponds to n s in the 0.94 — 0.96 
range). The amplitude of the vertex value, and therefore the coupling amplitude of 
test fields, is found to be amplified at the end of inflation by a factor that can be as 
large as 7 showing that quantitative results obtained in the de Sitter limit could be 
quite misleading! 




20. 30. 40. SO. 60. 

Ne 




0.8 I — £ _ 

20 30 40 50 60 

Ne 



Figure 2. Behavior of the effective vertex value as a function of the number of 
cfolds from horizon crossing. The value of the vertex is given in Log scale in units 
such that \/Hq = 1. 

Top panel: power law inflation. The solid line corresponds to v = 50.; the long 
dashed line to v = 41 and the short dashed line to v = 35. The dotted line to the 
de Sitter result. 

Bottom panel: inflation from a massive field. The solid line corresponds to a case 
where the total number of efolds is 50 whereas the long dashed line corresponds 
to N e = 60. The short dashed line corresponds to the de Sitter background case. 



4-3. The case of a massive field 

We pursue the exploration of the dependence of the results on the background 
dynamics with a more realistic, and data favored, inflationary model, that is inflation 
produced by a massive field. 
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In this paragraph we assume that the background evolution is derived from an 
infiaton potential Vi„f (<p) such, 



(4.14) 



where ip is the inflaton field. In order to simplify the calculation and to obtain closed 
forms, we assume that the slow roll solution can be used during the whole period of 
interest for our calculations. Simple numerical investigations show that this is indeed 
a good approximation. We then have, 



ip{t) = ip 



M P 



t 



(4.15) 



where t is the physical time, ipo is the value of <p at horizon crossing for the scales 
of interest and M P1 is the Planck mass. The time evolution of H can similarly be 
expressed as 

H(t) =H Q - Hft 

with 

1 



(4.16) 



#0 = 



m 



-ipo, and H 1 = ——. 



(4.17) 



The initial value ipo is determined by the number of foldings one imposes on the 
inflationary period. Within those approximations, we have, N e = ip^/AM^ and the 
time at the end of inflation is given by i on d = 2N b /Hq. 
We then obtain, 



9f (t,f)) 



1 



m 



exp 



Erfi 



■9HS 

2 TO 2 

H(t') 



Ho(t + t')-^$(t 2 +t' 2 ) 



Erfi 



~H(t) 
V2m 



(4.18) 



M;§,2;. 



(4.19) 



V2m 

where Erfi is the imaginary error function and the vertex value is finally given by 

H o 

where 2-F2 is the hypergeometric function. On the bottom panel of Fig. 2 we show how 
the effective vertex value evolves with time for two different hypothesis for N e . Once 
again, we see that the amplitude of the Non-Gaussian vertex value is large compared 
to the dc Sitter case and amplified by a factor slightly smaller than 10. 

What is then to be concluded for the ki dependence of the effective /nl (or tnl) 
parameters? This explicit dependence cannot be derived from the resolution of the 
classical system but requires the full resolution of the quantum system. However the 
results obtained here suggest that the wave-mode dependence of the non-Gaussian 
parameters is suppressed by a factor of about 7 compared to the de Sitter case. 
This implies that the shape dependence is at percent level whereas the global scale 
dependence is at sub-percent level. 
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5. The global structure of the correlation functions at tree order 

These calculations can obviously be pursued in a perturbative expansion, in particular 
for the calculation of higher order correlation functions. Unsurprisingly one expects 
a close relation between the quantum calculation and the classical calculations. More 
precisely we expect the leading order in g in the quantum calculations to be identical 
to the classical calculations. 

In the following we illustrate these calculations with the four-point function in 
case of a cubic potential. 



5.1. Quantum calculation 

Let us start with the expression of the leading term of four-point correlation function 
in case of a cubic potential. It naturally involves two vertices and it is therefore given 
by the expression of the expectation value of (xki • • • Xk 4 )c at second order in 
As an application of the relation (2.21) it is thus given by, 



(Xki ■ • ■ Xk 4 )c 



a* fa) 



Arf / dry" 



x <0| [[u kl ...u k4 ,tf( 7 >fa')l ,# (/) fa")l |0) (5.1) 



Replacing in this expression by its value leads to, 



(Xki • ••Xk 4 )c 



A 2 (2^) ; 



Sd X>i 



x fdr/atf) [ V dr/'afa")[G fel fa, v')G k2 fa,*/)- c.c] 

x [G k: Mv")GkAv,v")G\ kl+ ^\(v' ,V") - cc] + sym. (5.2) 

where c.c. stands for complexe conjugate and where sym. means that all symmetric 
terms should be included (there are 6 such contributions in total). When the G 
functions are replaced by their expression in terms of their real and imaginary parts, 
the result scales either in P 3 g 2 or in Pg A . Wc naturally expect the first type of 
terms to dominate for super-horizon scales. These are those that we consider in the 
following. That leads to the superhorizon expression of the correlation function, 

sh _ AW / 



(Xk! • • - Xk 4 ) 



a 4 (v) 



-So 




x fdf/afa') f dr/'afa") [P kl fa, rf) 9ka fa, r/) + g kl ( V , r,')Pk 2 ] 

x [Pk 3 (v,v')gk 4 (v,v")P\k 1+ k 2 \(v',r]") +sym. 

+PkMv")PkMv")g\ M2 \W,v") +sym.] . (5.3) 

In this expression the sym. terms account respectively for 11 and 5 other terms. As 
can be seen there are two types of contributions. Each one corresponds to a different 
term - and graphs - in the classical picture. This is what we show explicitly in the 
following. 
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Figure 3. Formal diagrammatic representation of the expansion (5.4). Each term 
has a specific order in the initial free field expression. This order is represented by 
the number of dots in the circle. The Wick theorem imposes that dots should be 

(i) 

associated by pairs. Each dot has value x\. : each circle represents the operator, 
Up/pl J d' j k<5rj(k — J^ki) where the explicit value of each vertex depends on the 
background dynamics as explained in the text. 




Figure 4. Diagrammatic representation of the four-point correlation function. 
The left side represents the standard quantum Feynmann diagram (note that in 
case of an expanding background, this representation is a oversimplified since 
there are different propagators to consider). The right hand side represents its 
super-Hubble limit where such a contribution is split into 2 parts. Those diagrams 
are obtained by gluing different terms of the expansion depicted in Fig. 3 following 
the Wick theorem (dots are associated by pairs). Once the integration over ki is 
done, each line then carries a power spectrum factor P%. (if). 



5.2. Classical calculation 

The calculation of (xki ■ ■ -Xk 4 )c in the classical picture is indeed given by a further 
expansion of x m terms of the initial Geld value, e.g. 

Xk = x k + x k + x k + • ■ ■ ( 5 - 4 ) 

where x is linear in the initial field, x is quadratic. The expansion is therefore 
different from the expansion of Eq. (3.11). As a result the expression of (%ki ■ ■ ■ Xk 4 )c 
at lowest order in A will be 

(x kl . . . xk 4 ) c = (xgxg^xSk + <xgx£°xgx£>« + s^ 5 - 5 ) 

where sym. accounts for respectively 5 and 3 other terms. The expressions of x can 
be expressed recursively in terms of the classical Green function and in terms of the 
free fields x ■ More precisely we have, 

xt\v) = -7^-r, f Vd 3 k' 9 k (v,v')a(v')u^( V ')u^,( V ') (5.6) 



and 



2(2tt) s 



A ' - f dr/d 3 k' g k {^ V <) a(r/) U « k ,(r/) 



2(2vr 

f , „ d 3 k ,; 

(2tt) 3 / 2 



d^T^^'^^'Oa^O^U-^'O^-k"^') (5-8) 
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It is straightforward to see that the two contributions in (5.5) precisely correspond to 
the 2 terms in (5.3). 

In the super-horizon regime, it is furthermore natural to take the super-horizon 
limit behavior for those quantities. That implies that 



V2 / 4 



(x kl ...X k ^ H = ^y(2 7 r) 3 <5 D 

Pk 1 Pk 3 P\ ]il +u 2 \ ^ dr]'g Q {r],r]')a 3 (»/)) + s y m - 

+P kl Pk 2 Pk 3 a( V ) I* dv'9o(v,v')a 2 (r)') f dr,"g ( V ', r/>V) + sym. (5.9) 

where the power spectrum are to be taken at equal time r\. Again one should properly 
take into account the symmetry factors. This expression further simplifies in, 

(Xk, • • • Xk 4 )f = ^ P kl Pk 3 P\^ 2 \ + sym. 

+ ^P kl P k2 P k3 + sym. (5.10) 
where u e s. is given by Eq. (4.1) and where z/3 is, 

^3 = ^r r dv'goiViV'Mv'fveS. (??')■ (5-11) 



0(77) 

The classical picture naturally leads to the construction of a whole set of such vertices 
that can be built recursively. This is what is shown below. 



5.3. The tree structure 

The general structure of the correlation functions of the classical problem taken at 
leading order is well known. Similar structures have been encountered in other aspect 
of cosmological instability growth (see Refs. [34, 35, 36]). It corresponds to a tree 
structure where each vertex point corresponds to a given order of \ in the expansion 
of Eq. (5.4). To be more precise, 

<x kl ...xOc= £ (x k p 1 l) ---x k p ; ) )c (5.12) 

{Pi},£iPi=2(n-l)) 

where the sum is to be done for all values of {pi} that are integers and with the 
constraints that ^ = 2(n — 1). This constraint comes from the fact that in order to 
connect n points, one needs at least n — 1 links and therefore the order of the product 
is at east 2(n — 1) whereas terms with more than n — 1 links are subdominant and are 
therefore not considered here. When represented in a diagrammatic expansion, as in 
Fig. 4, all contributing terms of Eq. (5.12) then form a tree. 

For the problem of interest here the calculation further simplifies since each order 
of x m the super-Hubble limit is given by a simple vertex value multiplied by a 
convolution of the zeroth order terms, 

^) = ^=ir/ fe ( k -E k «)^ 1 ...d>k l , 

xxS(v)---X^(v)- (5-13) 
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where v 2 — v e ft.- in other words the super-Hubble classical properties of the fields will 
be determined, at tree order, by a discrete set of numbers: the vertices that appear in 
all those constructions. This result extends the form given in Eq. (5.11) for the four- 
point function. All high order correlation functions form then a genuine tree structure 
characterized by the expression of the power spectrum and the vertex values. 

And because the vertices are all independent of k, it is actually relatively easy 
to derive them at once. Consider for that the related system defined by the motion 
equation, u = a(r))xo{v)> 

<(V) - ^jY«oM = -« 3 (/KWWll (5-14) 
aWJ d X 

which obviously exhibits the very same late time growing mode and same Green 

function*. We thus have, 

&\V) = ^ (5-15) 

so that the solution of this equation, taken at the time of interest, is nothing but 
the vertex generating function)!. This identification was already pointed out in [10] 
although with a looser connexion to the quantum field calculations. Solving (5.14) 
might prove difficult in general. This equation can however be solved numerically, or 
even analytically in some peculiar cases. It is particularly interesting to rewrite the 
motion equation (5.14) in time variable t and field variable xo, such that 

dV 

Xo + 3ff X o = -^[xo]. (5.16) 

The two equations are clearly equivalent. More precisely we define Q{t, \o) as the 
solution of (5.16) with the initial conditions Xo(*o = 0, r) = r and xo(^o = 0,r) = 
which corresponds to set Q to the growing mode of the evolution equation. The 
generating function Q (r) of the vertices is then nothing but 

Q(t) ee £ U\v)Y = Xo(W,r). (5.17) 

v 

where t cn d is the time at the end of inflation. 

We present in Figs. 5 the resulting shapes of the generating functions Q{f) for 
different cases. They are compared to approximate forms we derive below. 

The obtention of approximate forms is based on the observation that as long 
as the interaction term is small, the first term of the left hand side of (5.16) is 
negligible compared to the second (e.g. the time evolution if xo is weak). Within 
this approximation the evolution equation of xo is simply 

3#Xo - ^ [Xo] ■ (5.18) 
which can be formally solved in 

dxo At 

* At this stage we obtain what a direct application of the SN formalism would have given but with 

a very specific interpretation and as a mean to provide a controlled approximation of the general 

expression of the fourth cumulant given in Eq. (2.26) and of tree-order correlators in general. 

tt These generating functions can be used to build local PDFs (sec [36] for details on this construction 

in general and [10] for this peculiar case) but it is beyond the scope of this paper to present such a 

construction. 
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Figure 5. The shape of the vertex generating function S(t) for a de Sitter 
background (solid lines) compared to the approximate form of Eq. (5.21), dotted 
line, for a cubic potential (top) and quartic potential (bottom). The results 
are also compared to the power law inflation case (dashed lines). The curves 
correspond to \N e = 0.1 and N e = 60. The power law exponent corresponds to a 
spectral index of n s = 0.95 



For a potential in V(x) — A% p /p! we are left with 



G(r) 



1 + Xt p 



-2_P 
(P 



-2 r^ d dt 



(5.20) 



For a de Sitter background we then have 

p - 2 A N e 



Q dS (r) 



1 



rP- 2 



after i end 



(p- 1)! 3H 2 

to has been replaced by N e /H. For a power law inflation it leads to 
p - 2 A [exp(2JV e /i/) - 1] 



(5.21) 



1 



rP- 2 



(P-l)l 



(5.22) 



Interestingly, for this type of potential, the structures of the r dependence are 
the same and independent on the background dynamics. They differ only by the 
value of the lowest order vertex, v c s, . It means for instance that for a cubic potential 
v-i scales like v\, e.g. it grows like is the de Sitter limit. In other words, to a 
good approximation, the entire background dependence of the non-linear structure is 
encoded in the effective vertex expression that in the dependence given by the r.h.s of 
Eq. (4.4). 
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5.4- Conditions for large isocurvature fluctuations revisited 

In paragraphc 2.2 conditions were given on A for isocurvature perturbations to exist at 
horizon crossing. The result of the previous paragraph now suggests some validity rules 
for the perturbation theory calculations to be valid during the super-Hubble evolution. 
Indeed the previous calculations make sense only if of the nonlinear transform Q (r) is 
still close to a linear transform when its argument r is of the order of its typical value 
at horizon crossing, that is H. 

For a de Sitter background it implies that the dimensionless quantity that appears 
in the expression of (? ds (r), e.g. A iV e iJ p_4 /3, is less than unity. It simply states that 
v^H 2 should remain small. The condition is now more stringent than the one derived 
in paragraphe 2.2. For other backgrounds we can then exploit the fact the functions 
Q are very similar once expressed in terms of v^t 2 . This suggests that the validity 
conditions for other backgrounds simply reads v-^H 2 < 1 which introduces an extra 
factor 7 or so. The conditions is slightly more stringent than for de Sitter case, it now 
reads A < l/(7N e ) for a quartic potential, but leave significant room for the x to 
develop large fluctuations and yet develop significant non-Gaussianities. 

6. Conclusions and comments 

We have explored the dependence of the evolution of the high order correlation 
functions of a test scalar field with minimal coupling on the background dynamics. The 
exact resolution of this problem requires involved quantum field theory calculation and 
there are only a limited number of cases where the calculations can be fully completed. 
They correspond to a de Sitter background and to some specific cases of power law 
inflation. Some of these results were already known but a novel and effective method 
is presented here that allows to perform those calculations efficiently. The explicit k- 
dependence of the bispectrum and trispectrum is derived for such cases for respectively 
a cubic or quartic potential. 

It is furthermore explicitly shown how the super-horizon evolution is equivalent 
to the mode coupling evolution of a classical system. It is then possible to derive 
the super-horizon evolution of the amplitude of the higher order spectra for arbitrary 
background. The result is encapsulated in Eq. (4.4). This expression is computed for 
different classes of models such as general power law inflation or background of chaotic 
inflation. It is found that when the background evolution departs from a de Sitter 
behavior the coupling amplitude is amplified. For realistic cases, we found that this 
amplification is about a factor 7. Those results suggest that the evolution is actually 
minimal when the background is de Sitter. 

It is to be noticed however that the resolution of the related classical system does 
not permit, even within a perturbation scheme, to compute the exact ki dependence 
of the effective tnl and /nl parameters of the iso-curvature perturbations. Only in 
specific cases, where the quantum evolution can be computed, are those results known. 
If the background is close enough to de Sitter we found that the shape variations can 
induce variation of tnl or /nl parameters up to 7 percent when the number of efolds 
is about 50 from horizon crossing to the end of inflation. The scale dependence (e.g. 
dependence with respect to scale for a given configuration) is precisely given by 1/N e 
where N e is the number of efolds. Results obtained for a different background suggest 
that these dependence are about 7 times smaller for the inflationary backgrounds 
investigated in this paper. Those dependences are expected to be transferred to those 
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of the metric fluctuations for certain models of multiple- hybrid inflation. Note however 
that for models where the transfer of modes takes place sooner, that is after a much 
smaller number of cfolds, the fc-dependence can be made arbitrarily larger. Its precise 
value depends however on the details of the model and needs a full quantum calculation 
to be derived. 

We finally show that when one restricts the computation to tree order, once super- 
Hubble scales have been reached, the high-order correlation functions follow a genuine 
tree structure (see text for details) with /c-independent vertices. Correlation functions 
are then entirely encoded in the amplitude of the two-point correlation function and 
in the generating function of the vertices that can be explicitly computed for any 
backgrounds. Examples strongly suggest that, to a very good approximation, the 
background effects are entirely encoded in the expression of the lowest order vertex so 
that higher order correlation functions, when estimated for super-horizon limits, can 
indeed be built from this vertex value only. 
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